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We propose a scheme of searches for the Z' gauge boson as a virtual state in scattering processes at 
the Tevatron taking into account model-independent relations between the Z' couplings to fermions. 
We integrate the Drell-Yan process cross setion to construct two-parametric observables, which are 
suitable for Z' searches in the pp -> i+r process. The observables allow to constrain the Z' 
vector and axial-vector couplings to SM fermions in a general parameterization with non-universal 
Z' interactions with fermion generations. Also a one-parametric observable for searching for the 
popular leptophobic Z' boson is proposed. 
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I. INTRODUCTION 

A new heavy neutral vector boson (Z' boson) [l[ is 
a popular scenario of searching for physics beyond the 
standard model (SM) of elementary particles in modern 
collider experiments. Both the Tevatron and LHC col- 
I laborations try to catch the particle as a resonance in 
the Drell-Yan process. Observing no peak they conclude 
O | that the Z' mass is no less than approximately 1.79 TeV 
if one considers some predefined set of Z' models. 

Another approach is to search for Z' in processes where 
it manifests itself as a virtual state. This includes the 
processes with the so called low-energy neutral currents 
(LENC) mediated by a Z boson. The contribution from 
a virtual Z' state arises due to the Z — Z' mixing. For 
example, the data on parity violation in cesium can be 
used to constrain the Z' mass 0]. A combined analy- 
sis of data on atomic parity violation, inelastic neutrino 
scattering, and neutrino-electron scattering Q allowed 
to constrain Fermi-like couplings that effectively repre- 
sent Z'-mediated interactions at low energies J5|. Gen- 
eral review of low-energy constraints on the Z' boson is 
presented in Q and in Section 10 of Ref. 0- 

Significant amount of the Tevatron data is collected at 
the Z-boson peak at 66-116 GeV. At these energies the 
Z 1 boson also can manifest itself as an off-shell state, the 
Z coupling constants are influenced by the Z—Z' mixing, 
and these effects may allow to find Z' signals by fitting 
the experimental data. 

In order to select Z' off-shell hints, proper observables 
have to be introduced to amplify possible signal H . The 
signal generally means a deviation of some Z' parameter 
(i.e., a coupling constant) from zero at a specified con- 
fidence level. The more parameters interfere in the ob- 
servable, the weaker constraints on the parameters will 
be obtained. Thus, the key problem for off-shell Z' de- 
tection is to maximally reduce the number of the Z' cou- 
plings in the observable, which is used to fit the data. 



The ultimate scenario assumes a one-parametric observ- 
able. However, a two-parametric observable can be also 
useful and effective. For example, the strategy to con- 
struct observables driven by one or two parameters was 
successfully applied to analyze the final data of the LEP 
experiment leading to model-independent hints and con- 
straints on Z' couplings So, attempts of selecting 
possible Z 1 signals from Tevatron data seem to be per- 
spective. 

In this paper we investigate possibilities of constructing 
few-parametric observables for the Drell-Yan process tak- 
ing into account kinematics of the proton- antiproton col- 
lisions at v^S* = 1.96 TeV and model- independent para- 
meterization of the Z' couplings. Here we consider the 
case of a Z' boson with non-universal Z' couplings to 
fermion generations. The universality of couplings will 
be discussed in a separate paper, since it leads to a dif- 
ferent (reduced) initial set of coupling constants and, 
consequently, requires a separate procedure of construct- 
ing few-parametric observables. We conclude that two- 
parametric observables exist at energies corresponding to 
Z peak, and we obtain all of them. These observables can 
be used as a key to find possible signals of the off-shell 
Z' boson. Data fitting is a subject of a separate investi- 
gation and lies beyond the scope of the paper. 

The paper is organized as follows. In SectionHJwe pro- 
vide all necessary information on the low-energy Z' pa- 
rameterization for our calculations. Section IIIII contains 
specifics on Z' contribution to the Drell-Yan process, un- 
certainties, and kinematic variables suitable for hadron 
colliders. In Section HVl we construct the observables in a 
step-by-step manner. In Section [V] we briefly summarize 
and discuss the obtained results. The Appcndix n contains 
some supplemental numerical data. 



II. ABELIAN Z' COUPLINGS TO LEPTONS 
AND QUARKS 
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Being decoupled at energies of order of mz, the 
Abclian Z' boson interacts with the SM particles as an 
additional t/(l) gauge boson. Its couplings to the SM 
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fcrmions are usually parameterized by the effective La- 
grangian: 

C zff = \zjr [{vfi + J 5 a s fz l )cos8 
+(v f +7 5 a/)sin6»o] /, 



C z , fJ = -Z'J>f [(« / + 7 5 a/)cos0 o 



„5„SM 



(1) 



(Further details on the parameterization can be found 
in [13] ■) Here / is an arbitrary SM fermion state; a,f 
and Vf are the Z' couplings to the axial- vector and vec- 
tor fermion currents, respectively; 9q is the Z—Z' mixing 



angle; v 
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are the SM couplings of the Z-boson. 



This parameterization is suggested by a number of nat- 
ural conditions: 

• the Z' interactions of renormalizable types are to 
be dominant at low energies ~ friz- The non- 
renormalizable interactions generated at high ener- 
gies due to radiation corrections are suppressed by 
the inverse heavy mass l/mz> (or by other heav- 
ier scales 1/Aj <C 1/mz') and, therefore, at low 
energies can be neglected; 

• the Z' is the only neutral vector boson with the 
mass ~ mz> ■ 

At low energies the Z' couplings enter the cross section 
together with the inverse Z' mass, so it is convenient to 
introduce the dimcnsionlcss couplings 
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which are constrained by experiments. 

Below the Z' decoupling threshold the effective 17(1) 
symmetry is a trace of the renormalizability of an un- 
known complete model with the Z' boson, and it leads 
to additional relations between the Z' couplings [Tpj : 
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(3) 



where q u , qd, I, and v\ are an up- type and a down- type 
quark, a lepton, and a neutrino inside any fermion gen- 
eration, correspondingly, and a is a universal constant, 
which defines also the Z' coupling to the SM scalar fields 
and the Z—Z' mixing angle in (|T|): 
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As it was discussed in [T(|, the relations © cover a 
popular class of models based on the E 6 group (the so 
called LR, x-ip models). Thus, they describe correlations 
between Z' couplings for a wide set of models beyond 
the SM. That is the reason to call the relations model- 
independent ones. 



As a result, Z' couplings can be parameterized by seven 



independent constants a, v u , v c , v t , v e , v^ 



These 



parameters must be fitted in experiments. In a particular 
model, one has some specific values for them. In case 
when the model is unknown, these parameters remain 
potentially arbitrary numbers. 



III. ABELIAN Z' IN THE DRELL-YAN 
PROCESS 

At the Tevatron the most prominent signal of the 
Abelian Z' boson is expected in the pp — > scatter- 
ing process (Fig. [I}. The general idea of our approach 
is equally applicable both for dielectrons and dimuons in 
the final state. To be definite, we shall consider the di- 
electron case. Specifics concerning the dimuon final state 
will be addressed to in Section [TV] and in the Appendix^ 
The cross section of this process can be written in form 
of the partonic cross sections combined with the parton 
distribution functions (PDFs): 



d 3 



&AB 



dXqdXqdt 



= /^2 fq,A( X qiQ 2 )fq,B{Xq,Q 2 ) 



</■</ 



dt 



(5) 



where A, B mark the interacting hadrons (p or p) with 
the four-momenta kA, ks; f qt A{x q ,Q 2 ) is the PDF for 
the parton q in the hadron A with the momentum frac- 
tion x q (0 < x q < 1) at the factorization scale Q 2 . To 
access the parton distribution data, we use the MSTW 
The quantity a qq _ 



2008 package [11 



2 + e - is the parton- 
level cross section, which depends on the Mandelstam 
variable t = (p e + —p q ) 2 . All parton-level calculations are 
performed using FeynArts [l2| and FormCalc [l3| pack- 
ages. Hereafter, the hat over a variable denotes that this 
variable refers to the parton-level cross section. 
We define the PDF factor for each quark flavor: 

.fq,A( X qiQ 2 )fq,B{Xq,Q 2 ) = F q q(x q , Xq, Q). (6) 

The PDF factor and the parton-level cross section are cal- 
culated in the leading order (LO) in as, and <jab in the 
LO is obtained in this way. The ncxt-to-next-to-leading 
order (NNLO) corrections are then taken into account by 
multiplying uab by the NNLO K-factor, which is calcu- 
lated using the Vrap software 14 1 (see also ref. [H[). 
We also consider two kinds of uncertainties: 

• the PDF uncertainties A<j pbf . The MSTW 2008 
package provides 68% CL and 90% CL intervals. 
We consider the latter one; 

• the uncertainties due to the factorization scale vari- 
ation, Acq. To incorporate these uncertainties, 
we follow the common procedure and vary Q from 
\/I/2 to 2\/§, where s is a Mandelstam variable for 
the partonic level process: s = (p e + +p e -) 2 . 
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FIG. 1. Drell-Yan scattering process. A and B mark the 
interacting hadrons - p and p in the Tevatron case. 



The cross section then can be written as 

<7DY ± Af7pDF±A(TQ. 

The obtained triple-differential cross section provides 
full description for the Drell-Yan process. It is expressed 
in terms of three kinematic variables: x q , x q , and t. The 
shortcoming of these variables is that all three of them 
enter both the PDF multiplier and the parton-level cross 
section, since s is not an independent value (s = x q x q S). 

The quantities that are directly measured in exper- 
iments and used for event selection are the pseudora- 
pidities rj± and transverse momenta p^ of the final- 
state electrons. In the leading order in as the relation 
p^ = — p^ = px applies. The Mandelstam variables s, t 
and the momentum fractions expressed as 

M 2 



M 



4py cosh 



2 V+ - V- 



t 



M 



XV++V- )/2 



M 



i + e^+^^-y 2 ' 

('?++>?- )/2_ (7) 



Note, that x q and x q depend only on the sum of the 
electron pseudorapiditics, Y = (?] + + i]-)/2, while t is 
expressed in terms of the difference of the pscudorapidi- 
ties, y = (r/ + — ?y_)/2. The Y variable is the well-known 
intermediate-state rapidity, while y is related to the scat- 
tering angle in the qq — > e + e~ process and governs the 
parton-level kinematics (it can be found introduced in 
some textbooks, for example in [l6j). In this way the 
cross section is obtained as a function of M, Y, y: 



-J2F^M,Q 2 ,Y) d ^ 



dNWYdy 



dy 



s -(M,y). 



(8) 



Here, the Q-dcpcndcncc of F qq is shown just to indicate 
that we incorporate the uncertainties due to the scale 
variation in our analysis. 

Leading Z' contribution to the Drell-Yan process arises 
from interference between diagrams with 7* jZ and Z' 
intermediate states, resulting in corrections of order of 
0{g 2 ). The cross section reads as 



CTSM H 
-2 

a <Tn2 



av u cr av -T v u v e (r v v e 



Here a, Vf are the couplings defined in <j2j) , (j3]), and cr a 2, 
&av f j VvfVf, are the numerical factors that depend on M, 
Y , y. In this approximation there are six independent 
unknown quantities entering the Drell-Yan process cross 
section. In (J9j> the factors that include v u and v c arise 
only due to contributions of first and second generation 
fcrmions, respectively. The contribution from the third 
generation is neglected due to the nature of (anti)protons. 

Once again, we note that Y enters the PDF factors 
only, while y is included into the parton-level cross sec- 
tions only. This is a crucial point for our analysis, as it 
allows us to treat F qq and a qq ^ e + e - separately. There- 
fore, we can try to use any peculiarities in the M-, Y-, 
and y-dependence of the PDF factors and partonic cross 
sections to suppress some of the numerical factors in ^ . 
For example, in case after integration by one of the kine- 
matic variables over some specific region the factor er 0t , Ce 
appears to be much smaller than the other factors, we 
may neglect its contribution to the cross section and deal 
with five unknown parameters instead of six we had ini- 
tially. Of course, we assume that all the combinations of 
the Z' couplings in the cross section are of the same order 
of magnitude. The leptophobic Z' case, which seems to 
be a very popular parameterization nowadays, is treated 
separately in Sec. HVl 

In addition to the Z 1 couplings, there are another two 
unknown Z' parameters that affect ctdy- These are the 
Z' mass mz> and decay width Y z' ■ The latest data from 
the CMS and ATLAS indicates that Z' is heavier than 
1.79 TcV. This means, that for energies close to the Z 
peak the ctdy dependencies on mz* and T z> can be ne- 
glected, assuming that the Z' peak is far away. 

The Y and y values that we can investigate are limited 
by detector performance and conservation laws. From 
the condition < x q ^ q < 1 it is easy to obtain the in- 
dependent limits 



i n ^<r<i n ^. 

M ~ ~ M 



(10) 
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(9) 



The electromagnetic calorimeters of both Tevatron de- 
tectors, CDF and DO, cover the electron pscudorapidity 
range \r)±\ < 3.2 [13 . Therefore, 

|F|<3.2. (11) 

The limits for y are the same as for Y, 

This section can be briefly summarized by saying the 
following: the cross section of the Drell-Yan process con- 
tains six unknown linear-independent terms inspired by 
Z' boson. The cross section depends on three kinematic 
variables, which will be used in what follows to suppress 
some of the contributions from the unknown Z' parame- 
ters. This will allow us to amplify the signal of Z' that 
is possibly hidden in the Tevatron experimental data. 



IV. THE OBSERVABLE 

Of course, the most detailed description of a scatter- 
ing process is contained in the differential cross section. 
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But a possible Z' signal can be washed out by the in- 
terference between the six independent combinations of 
Z' couplings entering the cross section. In general, in- 
tegration by kinematic variables can leave this situation 
without changes. We need to pay special attention to the 
integration scheme to reduce the number of interfering 
parameters in order to make a successful data fit pos- 
sible. This scheme must derive benefits from kinematic 
properties of the cross section. 



A. Integrating by Y 

The intermediate-state rapidity Y enters the PDF fac- 
tors only. Let us study the M- and Y-dependence of 
F qq (M, Y) in Eq. dHJ. At any fixed kinematically al- 
lowed Y value F qq is a smooth monotonically decreasing 
function of M. Kinematic properties of F qq are different 
for each flavor but independent of Z' properties. So, the 
Y-dcpendcncc of the cross section can be utilized to sup- 
press the contributions of the second generation, i.e., the 
terms with av c and v c v e in Eq. 

We use the following integration scheme 



PDF Factor for it, GeV 



-l 



M=66 GeV 



^ dYW(M,Y)- d3<TUY 



l-Y m dYdMdy 
with a simple piecewise-constant weight function 
/ A(M), < \Y\ < Y u 



W(M,Y) 



Yi < |Y| < Y„ 



(12) 



(13) 



In Eq. (|12[) o\ denotes the value obtained by the inte- 
gration of the triple-differential cross section <tdy by Y. 
In fact, we integrate the PDF factor in Eq. ([8]): 



Fgg(M) =2 dM W{M, Y) 
Jo 

xF qq (M,Y), 



(14) 
(15) 



So, in in this subsection we will study F qq (M) for dif- 
ferent quark generations. An example of this integra- 
tion scheme is illustrated in Fig. O Note, that the Y- 
distribution for the Drell-Yan cross section is symmetric. 

The Y m value is some positive boundary chosen for 
the Y integration region. The exact maximal Y m could 
be determined from Eq. (fTU|) . However, this would be 
technically inconvenient. First, it depends on M. Sec- 
ond, it leads to difficulties in usage of the K-factor. We 
calculate the K-factor as 



K{M,Y) 



T NNLO 

7 SM 
LO ' 
°SM 



(16) 



For different considered regions of M values and for Y 
close to boundary values (|10p the numerical uncertainty 
of the calculations becomes large for Og^. Because of 



0.020 - 




FIG. 2. Plot illustrating the integration scheme with the 
weight function from Eq. (|13[> . The vertical dashed line rep- 
resents Y 1 . In this particular case Y 1 is set to 0.75, and the 
upper integration limit Y m is 2.35. 



this the K-factor becomes a non-monotonic fast-varying 
function of Y and, therefore, cannot be used to improve 
the new physics contributions to the cross section. To 
avoid these difficulties, we chose a somewhat lower value 
of Y m being also independent on M. 

First, let us consider the unit weight function (i.e., 
A(M) = 1). In this case there is no effect from Y\. The 
integration limit Y m is set to 1.65, which corresponds 
to a bin bound at the DO detector (see, for example, 
[3]). The plots of F q q(M) for u, d, c, and s quarks 
are presented in Fig. [3] (a). These plots indicate that 
for M > 240 GeV the factors for the second-generation 
quarks amount to less than 1% of those for the first gen- 
eration. This leads to the conclusion that for these values 
of M we can use the standard integration with A(M) = 1 
neglecting the contributions from the second generation, 
(Jav c and (Jv c vi ■ 

Now consider the M values at the Z-peak. Both CDF 
and DO collaborations define limits of this region to be 
symmetric with respect to the Z boson mass. These lim- 
its are often set to either 66 GeV < M < 116 GeV or 
71 GeV < M < 111 GeV plUl. In the present paper 
the former alternative is used. Actually, the choice of spe- 
cific lower and upper limits does not affect our results. 
There are only two general requirements: the limits have 
to be symmetric with respect to mz and large enough so 
that we could set all quark masses to zero. 

In Fig. [3]the plots for F qq {M 7 Y) versus Y at different 
M values are shown for u, d, c, and s quarks. The relative 
contributions of second generation quarks amount up to 
11% at M = 66 GeV and cannot be neglected. There is 
a qualitative difference between the PDF factors for the 
first and second generations. At some energies the factor 
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(a) 

FIG. 3. Plots illustrating the suppression of the contributions of second-generation quarks to the Drell-Yan cross section in 
different regions of M values and with different integration schemes. The plotted values are F U u, F dd , F C c, and F S s (not shown 
on the right plot. Because of the utilized integration scheme at some M is becomes negative, but its absolute value is even 
smaller than F C c) integrated by Y: a) The integration by Y is carried out over the region \Y\ < 1.65 with A(M) = 1; b) the 
integration by Y is carried out over the region \Y\ < 2.35 with A(M) from Fig. [5l where Y\ is set to 0.75. 



for u quarks is convex for Y close to zero (at somewhat 
lower energies this is also the case for d). This is due to 
the nature of a proton. 

For any given M value from the Z-peak region we can 
adjust the weight function in such a way that the factors 
F cS (M) and F sS (M) amount to less than 1% of each of 
the factors F uU (M) and F dd (M): 

F c5 , 3 - s (M)<0.01F uU>dS (M) (17) 

This is shown on Fig. [3] (b). Therefore, the contribu- 
tions of the second generation-fermions are suppressed, 
and again cr aCc and cr CcCc are excluded from ctdy- The 
weight coefficient A(M) can be determined for several M 
values and interpolated in the Z-peak region (see the Ap- 
pendix,). For our specific case A(M) is plotted in Fig. [S] 
Here Y m is set to 2.35, and Y\ is 0.75. 

As a result, we obtain the cross section a±, which de- 
pends on y and M and contains four linearly independent 



Z' terms instead of six: 

i-2 i — — 

o"i = (J i sm + a oi a 2 + av e a i aCe 

(18) 

Our next step is to use the remaining two kinematic 
variables, M and y, to get rid of another two unknown 
combinations of the Z' couplings. 



B. Integrating by M and y 

The difference of the pseudorapidities, y, enters the 
parton-level cross section of the Drell-Yan process, 
CT 9 q_j. e + e - , only and is irrelevant for the PDF analysis. 
The parton-level cross section depends also on M through 
four 'resonant' functions: 



MM) = 



(M 2 /m|-l) 2 + r|/r 



h(M) = 



(M 2 /m|-l) 



(M7m|-l)2+r|/m| 



(M*/m%, 1) 



(MVm|,-l) 2 + r|,/m|, 



h(M) 



M 2 TzTz'/(m 3 z m z ,) + (M 2 / 'm|)(M 2 / 'm\ - l)(M 2 /m|, - 1) 



p//m|-i) 2 + r|/m|; 



[(M 2 /m|, - l) a + r|,/m|, 



(19) 



Here mz,z> an d Tz,z' denote the masses and the widths case we do not care about the specific values of the Z' 
of the Z and Z' bosons. We investigate the energy region mass and decay widths. But at this point for numerical 
close to the Z boson peak. As it was noted earlier, in this calculations we arc going to set specific values for mz> 
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FIG. 4. Plots for F q q(M, Y) versus Y at different M values. The uncertainties that arise from the PDF errors and factorization 
scale variation are also shown (see Sec. IIIII) . 
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FIG. 5. Weight coefficient A(M) used for integration over the 
Z-peak region. Y m = 2.35, Yi = 0.75. 



FIG. 6. Plots for the resonant functions, which are given by 
Eqs. [H in the region 66 GeV < M < 116 GeV. 



and Tz'- Following the recent LHC results @], we set 
mz> to 1.8 TeV and assume the decay width to be 10% 
of the mass. That is, we use some asymptotics of //, and 
/ 3 at M < m z > ■ 

As it can be seen from Fig. [51 the /i function is dom- 
inant. The functions /2, f'i are odd- like with respect to 
AI = mz , and the function is small. As a consequence, 
after integrating by M over the discussed symmetric Z- 
peak region the functions f2, /;>, and fs arc negligible 
compared to f\. We are going to use the discussed fea- 
ture in what follows. 

When investigating the M -dependence of the hadronic 
cross section <j\ , we deal not with the resonant functions 
themselves, but with their products with the PDF fac- 
tors. The general form of o\ can be written as 

(7i - crj SM = [ a {M) tanh 2y + b(M)} , (20) 

cosh y 



where a(M) and b(M) are some functions that include 
the unknown couplings a, v u , and v e . The A/-dcpcndence 
arises from the 'resonant' functions multiplied by F q q(M) 
from Eq. (fl~4")) . From the plots in Fig. [3] we can con- 
clude that the factors F q g(M) are smooth, monotonic, 
and slowly-varying in the considered region. Therefore, 
we stress that all the discussed properties of fi, / 2 , 
and /3 are generally maintained, when these functions 
arc multiplied by F q q(M). 

Naturally, //> and do not enter the SM part ctism- 
There are four factors entering the Z' contribution: a x a 2 , 
ciac c , cias„, and cr iCuCe [see Eq. I|18p]. The factor 
cic„© c does not depend on fi, and, therefore, according 
to our discussion of properties of the 'resonant' functions 
we may eliminate it by the straightforward integration by 
M over the Z-peak region (66 GeV < M < 116 GeV). 
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FIG. 7. Plots for (a) the Z'-related factors and (b) 02 sm from Eq. (|21[) . The uncertainty bands are also shown. 



The resulting value is denoted cr 2 : 

C2 - cr 2 SM = y d-M (tJi — <Ti sm) 

cosh 2u . , „ , . 
= — - T ^(atanh2y + 6), 
cosh y 

02 = (T2SM + a 2 a 2 a 2 + aV e (J2av c + 0,V u <J2av u , 

a = J dMa(M), b = J dMb(M). (21) 

The factors cr 2 SM, 02a 2 , cr 2 as e , 02as u , and 02b„b c are 
plotted on Fig. It can be seen that cr 2 s u s c is negligibly 
small compared to the other three factors indicating that 
our assumption is relevant. 

We are not concerned about er 2 sm at the moment and 
shall turn to investigating the y-dcpendence of the Z'- 
rclated contribution presented in Eq. (p?Tj) . The behavior 
of the tTac„ factor is governed by its odd part, while the 
<7 5 2 and <Jav c factors are obviously dominated by their 
even parts. From the plots on Fig. [7J (a), one can con- 
clude that it is possible to suppress one of the three fac- 
tors by integrating the cross section by y over a symmet- 
ric region. Remember, that the integration limits for y 
are the same as for Y . In our case 



2.35 < y < 2.35. 



(22) 



For example, we can integrate them with a pieccwise- 
constant function 



w(y) 



x, y > 0, 
1, y<o. 



(23) 



Here x is some real number. The sign of x is chosen 
depending on which specific factor we want to suppress. 



The resulting observable a* would be a somewhat mod- 
ified forward-backward scattering asymmetry: 



dyu){y) <t 2 . 



(24) 



We propose an approach that is a bit more general. 
The weight function that we use has a structure similar 
to the one in Eq. ([2TJ]): 



ui(y) = tanh 2y + k. 



(25) 



Just like the Z' contribution to cr 2 , this is a sum of odd 
and even functions of y. Here k is a numerical constant. 
We will adjust its value so that the contribution of one 
of the remaining three factors becomes negligible when 
integrated by y. 

After the integration we obtain 



-'SM 



dy (tanh 2y + k) er 2 
tanh-i 



12a + bk 4- 



1 



cosh y 



-4a arctan(tanh y) 



a 



'SM 



(26) 



Note, that due to the symmetric integration region only 
the even part of the function w(y)cr 2 survives. The factors 



'SM' 



and er~ arc linear functions of k: 



ctsm = (3.40 + 63.5 k) pb ± (0.39 + 5.4 k) pb, 
cr f *2 = (0.354 - 12.2 k) nb ± (0.003 - 1.2 fc) nb, 
cr* 0u = (0.468 + 3.89 fc) nb ± (0.009 + 0.17 fc) nb, 
a* a - e = (7.12 + 0.802 fc) nb ± (0.52 + 0.068 fc) nb. (27) 
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TABLE I. Couplings entering each of the two considered ob- 
servables, together with the corresponding values of k, the 
SM contribution o"sm> an d the factors o"J 2 , viv u , and er| Se . 



couplings 


k 


CT SM>P b <2,nb 


°bb„ > nb <s e . nb 


a 2 , av u 


-9 


-569 ±48 111 ±10 


—34.5 ±1.5 suppressed 


a 2 , av e 


-0.12 


-4.23 ±0.26 1.82 ±0.14 


suppressed 7.02 ± 0.52 



Let us construct an observable that is suitable for fit- 
ting of the axial-vector coupling a and the coupling to 
the up-quark vector current, v u . That is, the factor 
has to be suppressed. We choose the suppression criteria 

KU< 0.011^1, K.|<0.01|o* j J (28) 

to calculate k in Eq. (|27|). Overlap of the intervals ob- 
tained from the upper and lower bounds for factors gives 
the resulting interval —9.18 < k < —8.55. If we set fc = 
-9 in Eq. (|27|) . the resulting observable will contain only 
two unknown Z' parameters: 

a * = °sm + a 2 <J* 2 + av u cr^, u , 
a* SM = -569 ± 48 pb, 
cr? 2 = 111 ± 10 nb, 
cr*^ = -34.5 ± 1.5 nb. (29) 

This specific observable allows us to perform fitting of 
the a and v u couplings. 

There are two other possible observables in this ap- 
proach: the one with suppressed cr| Cu and the one with 
suppressed <j* 2 . However, the latter case cannot be real- 
ized in our scheme with suppression factor 0.01, because 
the intervals obtained for the lower and upper bounds 
from ([27]) do not overlap. Therefore, one has either to re- 
quire weaker suppression in Eq. (|28[) or to narrow down 
the margin of error reducing the confidence level. Fur- 
thermore, this observable contains three Z' couplings as 
opposed to two couplings in the case when <r a 2 or Cg C 
is suppressed. The mentioned flaws make this observ- 
able less attractive for data fitting, and we refrain from 
discussing it in the rest of our paper. 

In Table [J we present the combinations of couplings 
that enter each of the proposed observables, together 
with the corresponding values of k and a* 2 , a~ , and 
tTg C . Note, that we choose certain k values from the 
corresponding intervals. 

The model-independent analysis of the LEP II 
data [Io| resulted in obtaining upper bounds for a 2 and v 2 
at 95% CL, both of order of 10" 4 . From Figs. [3(a), (b) 
and Table U (see also [2l|) it can be seen that these upper 
bounds are too large, since when substituted into Eq. [5T] 
they lead to a large deviation from the SM, which is not 
confirmed by any of the experimental data. Therefore, 
we may expect at least some significant improvement of 
the LEP-motivated bounds. 

Neither LEP data nor Tevatron or LHC data shows 
any explicit indications of the Abelian Z' . This provides 



motivation to investigate models with the so called lep- 
tophobic Z' [22||. In these models Z' boson couplings 
to the SM leptons are strongly suppressed compared to 
the quark couplings. Among other things, this param- 
eterization allows to explain deviations of the precision 
electroweak data from the SM by introducing Z' with the 
mass close to mz [H[. From the Lagrangian in Eq. ([T]) 
and the relations in Eq. ([3]) it follows that in the lepto- 
phobic case vi, a; , and a q are small compared to v q , and 
the leading Z' contributions to the cross section are 

+av c <Jav c + v c v e av c % + 0(a 2 ,av e ). (30) 

After applying all the integrations discussed in Section 
ITVl we end up with the observable where only the term 
av u a^ survives. This observable is one-parametric: 

CT * = cr SM + ««« cr ac„- (31) 

The numerical values are the same as in the second line 
of Table [J 

Our results obtained for the diclcctron case can be 
easily recalculated for dimuons, taking into account the 
difference between detector pscudorapidity coverages for 
electrons and muons. For example, the CDF Collabo- 
ration detects electrons with maximum pseudorapidity 
he | = 3.2 H3 (the fiducial region is \r) e \ < 2.8). For 
muons this value is |?7 M | = 1.5 [24| . therefore, the value 
of Y m for the pp — > process is lower than for the 

dielectron case. This leads to different weight functions 
and k values, which are presented in the Appendix, 

V. DISCUSSION 

The data analysis performed by the LHC and Tevatron 
collaborations resulted in setting model-dependent lower 
bounds on the Z' mass. In that analysis only the high- 
energy region of the Drell-Yan cross section was consid- 
ered. In our paper we present a different approach that 
allows to search for a Z' signal in the pp — > l + l~ pro- 
cess at the energies near mz- In this region the most 
important contributions at the Z peak come from the 
Z — Z' mixing angle and Z'-induced contact couplings. 
The Z — Z' interference has to be taken into account, 
since it affects resonance shape as it was discussed in 
[25| . The approach utilizes the model-independent re- 
lations between the effective Z' couplings. Therefore, in 
case no signal is observed one would still be able to derive 
constraints for different Z' models and compare them to 
the ones presented in . 

The proposed prescription includes the following steps: 
1. The triple-differential cross section of the Drell-Yan 
process is expressed in terms of three kinematic variables: 
the mass of an intermediate state, M, the intermediate- 
state rapidity Y, and the variable that describes the par- 
ton scattering subprocess, y. This cross section contains 
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six unknown combinations of the Z' couplings: a 2 , av u , 
avi, av c , v u v h and v c Vf, 

2. The cross section is integrated by Y over the 
symmetric region [— Y m ]Y m ] with the weight function 
W(M,Y) defined in Eq. (HQ. The integration limits 
have to be determined for each specific final state (e + e~ 
or /i + p,~) and detector individually. The function A(M) 
has to be adjusted in such way, that the PDF factors 
for the second-generation quarks, F cS (M) and F sS (M), 
amount to less than 1% of the PDF factor F U u(M). As 
a result we exclude av c and v c vi from the cross section; 

3. Integrate the cross section by M over the Z bo- 
son peak region: either 66 GeV < M < 116 GeV or 
71 GeV < M < 111 GeV, or any other region with 
bounds symmetric with respect to mz- These bounds 
have to be large enough, so that one could neglect the 
masses of the u, d, c, and s quarks compared to M. 
This integration suppresses the contribution of v u vi to 
the cross section; 

4. The integration by y with the properly adjusted 
weight function uj(y) from Eq. (|25|) allows to suppress 
either av u or avi. 

The obtained two alternative observables can be used 
in fitting the experimental data on the pp -> scatter- 
ing collected by the Tevatron collaborations. This allows 
to constrain the Z' vector axial- vector couplings to SM 
fcrmions. 

In case of the leptophobic Z' boson, there is a one- 
parametric observable containing the combination of cou- 
plings dv u . 

There is a large amount of data on leptonic scattering 
processes collected in the LEP and LEP II experiments. 
The second observable in Table U contains the coupling 
combinations a 2 and av e that also enter lepton scattering 
processes. It seems to be attractive for combined fits of 
the LEP and Tevatron data. 



Appendix: Numerical data for the e e and \x p, 
cases 

In this section we provide numerical values of the limits 
of integration by Y and y, weight functions W(M, Y) and 
uj(y), and the observable a* . 

Generally, our approach is applicable for any final 
dileptonic state. However, the detector coverage is differ- 
ent for electrons and muons. The detector performance 
affects the proposed limits of integration by M, Y, and y. 
While the considered range for M is well-covered by both 
CDF and DO detectors, the detector-imposed limitations 
on the Y and y variables need a closer look. 

Both the DO and CDF Collaborations detect electrons 
with maximum pseudorapidity \r] e \ = 3.2 (l7| . There- 
fore, the proposed integration limits \Y\ < 2.35 and 
\y\ < 2.35 are applicable in case of the pp — > e + e~ pro- 
cess. However, the pseudorapidity coverages for muons 
are |^| < 2.0 and |r^| < 1.5 for DO [H and CDF [H, 
respectively. Because of this we have to appropriately 



adjust the limits of integration by Y and y. 

It is convenient to set Y\ = 0.75 for all three cases, 
dielectrons at DO and CDF, dimuons at DO, and dimuons 
at CDF. We choose the Y m value as shown in Table UH 
since it represents the integration limit both for Y and y. 
In Table HTT1 we present the values of A(M), which enters 
the weight function W(M, Y). 

TABLE 11. Numerical values of Y m for different final-state 
dileptons at DO and CDF. 





e + e 




> 






D0 


CDF 


Y 


2.35 


2.0 


1.5 



TABLE III. Numerical values of A(M) for different final-state 
dileptons at DO and CDF. 



M, GeV A(M), e+e" A(M), fi + pT 







DO 


CDF 


66 


-0.990053 


-0.929972 


-0.705828 


71 


-0.941569 


-0.884751 


-0.684274 


76 


-0.896807 


-0.844891 


-0.664090 


81 


-0.855355 


-0.809493 


-0.645150 


86 


-0.816857 


-0.777846 


-0.627341 


91 


-0.781010 


-0.749385 


-0.610566 


96 


-0.747549 


-0.723652 


-0.594736 


101 


-0.716242 


-0.700272 


-0.579775 


106 


-0.686890 


-0.678936 


-0.565613 


111 


-0.659313 


-0.659389 


-0.552186 


116 


-0.633356 


-0.641413 


-0.539441 



For interpolation we express A(M) as 



The values of oi, a 2 , and 03 are presented in Table ITVl 



TABLE IV. Numerical values of 01, 02, and 0,3 for different 
final-state dileptons at DO and CDF. 





e + e 


A 1 4 


> 






DO 


CDF 




0.225 


-0.212 


-0.049 


02 


146. 


53.4 


96.8 


CL3 


54.3 


8.35 


81.3 



The weight function w(y) from Eq. (|2"5"|) . as well as the 
rest of the prescription, is used for both dimuonic cases. 
For the DO case the resulting fc-dependent factors are: 

cr£ M = (2.84 + 50.8 k) pb ± (0.29 + 3.9 k) pb, 
(7^2 = (0.251 - 9.86 k) nb ± (0.002 + 0.85 k) nb, 
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cr*- u = (0.349 + 3.00 k) nb ± (0.004 + 0.11 k) nb, 
cr*^ = (5.54 + 0.648 k) nb ± (0.37 + 0.049 k) nb. 

For the CDF case we have 

crg M = (1.33 + 23.3 k) pb ± (0.12 + 1.6 k) pb, 
cr*2 = (0.090 - 4.57 k) nb ± (0.003 + 0.35 k) nb, 
a* Mu = (0.139 + 1.31 k) nb ± (0.001 - 0.03 k) nb, 
cr! Cfi = (2.36 + 0.302 fc)nb± (0.13 + 0.019 k) nb. 

Values of k obtained from suppression criteria (|28| for 



DO and CDF cases are shown in Table fVl 



TABLE V. Couplings entering each of the two considered ob- 
servables, together with the corresponding values of k, the 
SM contribution ag M , and the factors er* 2 , civ uJ an d cr^ ije . 



couplings 


k 


"SM> P b °fi2 ! nb 


crf_ , nb 


a-- , nb 








DO case (Y m = 2.0): 






a 2 , 


av u 


-8.5 


-429 ± 33 84.0 ±7.2 


-25.2 ± 0.9 


suppressed 


a 2 , 


av,j. 


-0.116 


-3.06 ±0.17 1.40 ±0.10 


suppressed 


5.47 ± 0.36 


CDF case (Y m = 1.5): 


a 2 , 


av u 


-7.8 


-180 ± 12 35.8 ±2.8 


-10.1 ± 0.2 


suppressed 


a 2 , 




-0.106 


-1.14 ± 0.05 0.575 ± 0.034 


suppressed 


2.33 ± 0.13 
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